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ABSTRACT 

vThe  prediction  of  torsional  buckling  of  stiffeners  is  a  necessary 
design  skill  in  the  planning  of  any  framed  structure  or  vehicle.  The 
likelihood  of  torsional  failure  is  increased  when  the  stiffener  is 
simultaneously  loaded  from  perpendicular  directions,  a  common  loading 
condition  for  the  frames  of  a  ship  at  sea.  This  thesis  examines  first 
the  case  of  a  single  lateral  load  on  I  beam  and  flat  bar  stiffeners 
and  then  the  instance  of  combined  lateral  and  axial  loadings.  An 
energy  method  is  usr  d  and  formulated  are  expressions  for  the  critical 
buckling  load  for  a  single  laterally  applied  force  and  a  formula 


relating  lateral  and  axial  leads  and  giving  an  indication  of  the 
relative  weakening  of  a  stiffener  caused  by  simultaneous  application 
of  both  forces.^ 

\ 
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NOTATIONS 


Length  of  stiffener 
Breadth  of  I  beam  flange 

Angle  of  rotation  of  stiffener  in  the  YZ  plane 
Amplitude  of  cosine  function  for  B 
Height  of  stiffener  shear  center  above  toe 
Warping  Constant 

Additional  lowering  of  the  load  P  due  to  the  load 
being  applied  a  distance  h/2  above  center 
of  the  web 

Young's  Modulus  =  2.068  x  10®  KN/m 
Shear  Modulus  =  E/2.6 
Height  of  stiffener 

Moment  of  inertia  of  stiffener  about  web  plane 
Local  rectangular  coordinates  at  any  point  x 
on  a  stiffener 

St.  Venant's  torsion  constant 
Dimensionless  parameter 
Moment  about  the  l  axis 
Lateral  load 

Critical  lateral  buckling  load 
Axial  end  load 


Euler  buckling  load 


Flange  thickness 
Web  thickness 

Displacements  in  the  x,y,  and  z  directions 

Strain  energy 

Work 

Rectangular  coordinates  originating  at 
centroid  of  stiffener 
Numerical  factor 
Membrane  strain 


Poisson's  Ratio 
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INTRODUCTION 


As  is  the  case  with  most  large  engineering  structures,  a  ship's 
hull  and  its  interior  dividing  bulkheads  and  decks  are  basically 
framework  patterns  of  stiffeners  covered  by  an  outer  shell.  It  is 
principally  this  lattice  of  stiffeners  which  determines  the  sti  *■  yth 
and  structural  stability  of  a  vessel.  Structural  soundness  and 
integrity  are  particularly  critical  for  a  ship  because  this 
engineering  structure  must  operate  in  what  can  be  the  harshest  natural 
environment  on  earth.  Very  seldom  are  the  frames  of  a  ship's  hull 
subject  to  loading  in  a  single  plane;  rather,  it  is  most  often  a 
combination  of  loads  caused  by  the  interaction  of  the  various  natural 
and  man-made  forces  from  different  directions  which  determine  the 
total  stress  on  a  member.  The  subject  of  this  thesis  is  to  develop  an 
expression  relating  the  interaction  of  simultaneous  lateral  and  axial 
loading  of  a  stiffener. 

Two  of  the  more  commonly  found  shipboard  stiffeners,  the  I  beam 
and  flat  bar  type,  will  be  examined  first  under  the  action  of  a 
centrally  applied  lateral  load  and  then  under  the  combined  actions  of 
the  central  lateral  load  and  a  per pindicularly  applied  end  axial  load. 
Concentrated  point  loads  were  chosen  because  they  represent  the  most 
limiting  and  potentially  dangerous  cases.  An  energy  approach  will  be 
used  in  all  instances.  First  the  strain  energy  of  the  system  will  be 
found  and  then  the  work  done  by  the  loading  forces.  By  conservation 
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of  energy,  the  variation  of  the  strain  energy  will  be  equated  to  the 
work  done.  From  the  resulting  expression  will  be  developed,  in  the 
case  of  the  single  lateral  load,  a  formula  for  the  critical  buckling 
load  and,  in  the  case  of  combined  loading,  a  formula  relating  the 
interaction  of  the  lateral  and  axial  forces. 

It  is  hoped  that  study  into  the  problem  of  stiffener  failure 
under  more  realistic  loading  conditions  such  as  this  may  aid  in  the 
formulation  and  upgrading  of  design  codes  used  in  civilian  and 
military  ship  design,  and  any  efforts  that  could  make  a  ship  safer  or 
lighter  would  certainly  be  worthwhile.  To  this  end  it  is  hoped  that 
the  ultimate  contribution  of  this  thesis  is  a  reasonably  accurate 
method  of  predicting,  given  a  value  of  axial  load,  how  much  lateral 
load  may  be  sustained  by  a  stiffener  without  buckling,  or  vice  versa. 


CHAPTER  1 


DESCRIPTION  OF  MODELS 

The  two  types  of  strengthening  members  that  will  be  considered 
are  the  I  beam  type  and  flat  bar  type,  to  be  known  respectively  as 
Model  I  and  Model  II.  The  coordinate  system  to  be  employed  with  both 
models  is  shown  in  figure  (1-1). 


Z 


Figure  (1-1)  Coordinate  System 


The  origin  of  the  coordinate  system  will  be  taken  at  the  center 
of  the  beam,  with  u,v,  and  w  being  the  displacements  along  the  x,  y, 
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1.1  Model  I 


Model  I  as  stated  has  been  chosen  to  be  a  simply  supported  I 
beam.  The  case  of  a  beam  with  symmetric  flanges  will  be  considered, 
and  a  lateral  cross  section  of  Model  I  is  shown  in  figure  (1-3). 


Figure  (1-3)  Model  I  Lateral  Cross  Section 


The  properties  of  the  Model  I  cross  section  are  listed  below: 


c  -  Height  of  shear  center  above  toe  (Flange 


thickness  t,f  considered  negligible  compared 
to  web  height  h) 


c  = 


h 


7 


Cw  -  Warping  constant 


Cw 


b 


3 


J  -  St.  Venant's  torsion  constant 


J  =  — ~  ( 2bt  f  3  +  htw3) 


Iz  -  Moment  of  inertia  about  the  web  plane 
(contribution  of  top  and  bottom 
flanges  only) 


Iz 


b3t . 


2  ( 
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) 


15 


»*  *.  *  V  >  »  •  »  '  ,  *  „*•»  .*»  A\N*  \  V* 


The  properties  of  the  flat  bar  cross  section  are  as  follows 


c  -  Height  of  shear  center  above  toe 


h 

C  ~  T 


Cw  -  Warping  Constant 


Cw  =  0 


J  -  St.  Venant's  torsion  constant 


I„  -  Moment  of  inertia  about  the  web  plane 


CHAPTER  2 


Model  I  -  TORSIONAL  BUCKLING  UNDER  POINT  LATERAL  LOADING 


2.1  Development  of  the  Strain  Energy  Equation 


The  case  of  a  simply  supported  I  beam  of  length  a  subjected  to  a 
point  lateral  load  P  in  the  center  is  illustrated  in  figure  (2-1) 


Figure  (2-1)  I  Beam  Under  Point  Lateral  Load 


The  general  equation  for  the  strain  energy  of  a  beam  of  length  a  under 
bending  and  torsion  is  given  by  (refer nee  1) : 


U  =  C  [""eI  (  d  %-)  2  +  GJ  (4^-) 2  +  E  Cw  (  I  )  2  I  dx 


Considering  only  that  portion  of  the  beam  to  the  right  of  an  arbitrary 
point  x,  figure  (2-2),  the  forces  acting  will  be  a  single  one  P/2  in 
magnitude  acting  at  x  =  a/2.  If  a  local  coordinate  system  with  axes 
j,  k,  £  is  assigned  to  the  point  x,  the  force  P/2  will  produce  a 
moment  about  the  &  axis  equal  to: 


P 

M  H  =  (  a/2  -  x ) B  (2-1) 


with  force  P/2  acting  to  the  right  a  distance  a/2  -  x 

Where  B  again  is  the  angle  of  twist  indicated  in  figure  (1-2).  It 
should  be  noted  that  B  is  assumed  to  be  sufficiently  small  so  that 
small  angle  approximations  for  sine  and  cosine  are  applicable. 


The  expression  for 


d2v 

dx^ 


then  is  (reference  1) : 


d2v 

dx2 


M, 


El 
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.  .  ^x 

For  B  the  simple  choice  B  =  Bo  cos  ~  will  be  made.  The  cosine 
expression  is  selected  because  it  is  assumed  that  the  stiffener  will 
be  mounted  such  that  it  is  restrained  from  twisting  at  the  ends  x  = 
a/2  and  x  =  -a/2.  This  is  a  reasonable  assumption  in  the  case  of 

stiffeners  fixed  between  decks  or  frames  of  a  ship. 


Utilyzing  the  appropriate  terms  for  M  ^  and  3, 


d2v  _  P  ,  a 

.  2  2EI  '  2 

dx  z 


x)  Bo  cos 


TTX 

a 


becomes : 


(2-2) 


Inserting  this  expression  into  the  formula  for  strain  energy  U  yields: 


U  = 


//2  PpW 

Jo  L4EI* 


2  2 

P  Bo  _ 2  iTx  ,  a  v2  ,  2  2  tx 

4ett  cos  t  1  r  +  GJ<— 1  sin  — 


+  ECw  Bo2  (  -  ) ^  cos2 

3, 


€ 


dx 


Carrying  out  the  integration  it  becomes: 

U  =  -pfsal  a3(  S-*  /-  )  +  GJ  (  -Sfl.)2  |  +  ECvj  BO2 


El 


192- 


2.2  Development  of  the  Virtual  Work  Equation 


a 

4 


Work  is  defined  as  the  movement  of  a  force  through  a  distance. 


In  the  case  of  Model  I  then  the  virtual  work  done  is  obtained  by 
multiplying  the  force  P  times  the  summation  of  all  the  components  of 


motion  in  the  Z  direction  for  elements  of  the  beam  between  x  =  0  and  x 
=  a/2.  It  is  not  necessary  to  sum  vertical  components  over  the  entire 
length  of  the  beam,  x  =  -a/2  to  x  =  a/2  ,  because  each  half  moves  the 
same  total  distance.  The  expression  for  the  vertical  distance  moved 
by  the  centroid  of  a  cross  section  of  the  beam  at  the  arbitrary  point 
X  is  ( reference  1 )  : 


B 


and  the  work  equation  then  becomes: 

a/2 

W  =  P  C  B  —  (  |  -  x  )  dx 

Jo  dX 

This  equation  would  suffice  alone  if  the  force  P  were  applied  at  the 
centroid  of  the  cross  section  at  X  =  0;  but  since  we  are  considering 
the  more  realistic  case  where  the  load  is  applied  at  the  outside 
surface  of  the  beam,  a  distance  h/2  above  the  centroid,  the  expression 
must  be  slightly  ammended.  Due  to  the  twist  of  the  section  at  X=0, 
the  load  will  drop  through  an  additional  distance  d  =  h/2(l  -  cos  B*.0  ) 
as  shown  in  figure  (2-3). 
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Figure  (2-3)  Additional  lowering  of  P  due  to  load 
being  applied  a  distance  h/2  above  centroid 

By  the  half  angle  relation, 

2 

1  -  cos  Bxeo=  2  sin  B,^/2 

=  2  (B»m/2  )  2 

—  B  k-o /2 

Therefore  d  becomes  approximately: 

d  -  (h/2)  (bL/2) 

=  (h/4)  Bo2  cos20 
=  (h/4)  Bo2 


This  terra  is  now  multiplied  by  P  and  added  to  the  work  equation  to 

give  a  more  accurate  expression  for  the  work  done. 
a/2 


W  =  P 


5. 


B 


d2v 


dx‘ 


(  2 — x)dx  + 


Ph  2 
—Bo 


(2-3] 


a/2 


=  P 


Bo  cos 


TTX 


2 

,  d  v  w  a  ,  ,  ,  Ph  „  2 
( - 2  ^  (~ ~x)dx  +  - Bo 


dx 


Again 


d2v 


2EI 


(~2~  ~  x  )  B  •  an(^ 


dx  z 

when  applied  the  work  equation  becomes: 


w  = 


2  2 
P  Bo 


a/2 


2EI 


i 


COS 


rx  (  a_ 

a  v  2 


-  x  )2  dx  +  7^-  BO2 


Integrating  this  expression  gives 


„  P2Bo2  3  ,  6  +  -2  v  ,  Ph  „ _ 2 

W  =  -  a  (  - —  )  +  -7-  Bo 


El 


96  ~ 


2.3  Determination  of  the  Critical  Buckling  Load 


To  determine  the  value  of  P  that  will  cause  failure  of  the 


stiffener,  the  critical  buckling  load  Per,  the  principle  of 
conservation  of  energy  will  be  employed.  This  concept  states  that 


when  a  member  is  gradually  loaded,  the  kinetic  energy  is  zero  and  the 
work  done  by  external  forces  is  equal  to  the  strain  energy  (reference 
2)  : 

6W  =  5U 


To  determine  6W  and  <$U  calculus  of  variations  with  respect  to  Bo  is 
applied  to  the  work  and  strain  energy  equations.  First  to  the  work 
equation : 

2 


W(Bo  +  6 Bo)  =  (Bo  +  ^Bo) 


2  3  2 

Pa  .  6  +  r*  Ph 

El  q  ^  2  4 

z  96 


1“  2  3 

=  W (Bo2 )  +  2Bo ‘ Bo  ^ 

2  z  96  ~  2 


2)  + 


oW  =  2 Bo  *  Bo  |  a——  (  6  ~V  )  +  ^-1 


fpi 

El. 


96" 


U 


And  in  like  manner  to  the  strain  energy  equation: 


3  2  - 1 

U  ( Bo+  ;  Bo)  =  ( Bo+  ''Bo)  2  |  2_|.  (_§±—)+  GJ(-)2  |  +  ECW(^-)2  |  j 


P 
0 


z  192" 


9  I  P2.3  £  +  ,2  -  2 

=  U  (Bo)  +Bo s  Bo  |  —  <- - j)  +  4a  (GJ  +  ECw( 


z  192- 
2  3  2  _2 

<5U=2Bo^Bo  |  |t^-(6  +  :  5-)+  7—  (GJ  +  ECw  ( 
EIZ  192” 2  4a 


tf 


^'2] 


The  two  variational  expressions  may  now  be  equated,  and  the  formula 
for  the  critical  buckling  load  derived. 


2  Boo  Bo  |  2£|_a_  (  2±I_)  +  _P££il 


] 


El 


96’ 


0 


2bo6bo  I  EfE^ai  (  5±i  )  + 
z  192^ 


„2  ,  ,  21 
4a  (GJ  +  EC"<i>  J 


2  3  2 

Per  a  ,  6+n  .  Per 

*  n  +  /r 


EI  192tt‘ 

z 


u  2  r 

•  7T 

u  n 

4a 


GJ  +  ECw  (— 


T>2]  *  • 


This  quadratic  may  be  solved  for  Per,  and  the  result  expressed  in  the 


f  orm : 


£ 


6+TT 


Per  -1  -  (  214.,  6  \/-Elz 


GJ 


+  1 


,  4 8  r 2  .2  h.  2  FIz 
'  2  j  V  GJ 


64- TT 


+  (■ 


192- 


)  (  1  + 


ECw 


6  + 


GJ 


<  i  >  > 


1/2  )  /EI  GJ 


(2-4) 


This  formula  corresponds  directly  with  equation  (6-18)  of  reference  1 
for  similar  geometry  and  loading  which,  in  the  notation  of  this 
thesis,  is: 

V  Vei 

z 


Per  = 


y Vei  gj 


The  expression  within  braces  in  equation  (2-4)  equates  to  the  non- 
dimensional  factor  y  of  (2-5).  Figure  (2-4)  shows  the  close  agreement 
between  the  values  of  Y  calculated  from  equation  (2-4)  with  those 
listed  in  tabular  form  in  Table  6-5  of  reference  1. 


CHAPTER  3 


MODEL  II  -  TORSIONAL  BUCKLIN3  UNDER  POINT  LATERAL  LOADING 


3.1  Development  of  the  Strain  Energy  Sguation 


Model  II  again  was  the  flat  bar  type  stiffener  illustrated  in 
figure  (1-4).  In  this  chapter  the  same  type  of  loading  situation  as 
that  of  chapter  2,  a  point  load  coincident  with  the  Z  axis  and  applied 
at  the  upper  surface  of  the  stiffener,  will  be  examined.  A  pictorial 
representation  of  the  loading  will  appear  identical  to  that  shown  in 
figure  (2-1) . 

The  equation  for  the  strain  energy  of  a  flat  bar  in  lateral 
bending  and  torsion  can  be  obtained  in  exactly  the  same  manner  as  that 
for  the  I  beam.  The  only  difference  is  that  the  warping  rigiditv 
term,  Cw,  vanishes.  This  is  because  the  equation  for  Cw  of  an  I  beam 
is : 

tfh2b3 

Cw  =  - 
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and  the  flat  bar  stiffener  can  be  likened  to  an  I  beam  with  the  width 
b  of  the  top  and  bottom  flanges  equal  to  zero.  The  strain  energy 


expression  then  is: 


a/2 


EIz  <^>2  +  GJ  ^ 
_  dx 


The  simple  choice  for  B,  Bo  cos  - —  ,  is  again  made  and  the 

,2  a 

expression  for  js  then  identical  to  (2-1)  : 

dx2 
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The  complete  strain  energy  equation  is: 
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which  when  integrated  becomes 
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3.2  Development  of  the  Virtual  Work  Equation 


The  work  equation  is  found  again  by  defining  the  work  done  as  the 
product  of  the  load  P  times  the  summation  of  all  vertical  components 
of  motion  over  one  half  of  the  bar.  The  expression  for  the  vertical 
distance  moved  by  the  centroid  at  any  point  along  the  X  axis  is 
identical  to  the  I  beam  case: 
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,2 

d  v 
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x) 


Also,  a  correction  term  must  once  more  be  added  to  the  work  equation 

to  account  for  the  load  being  applied  at  the  upper  surface  of  the 

plate  rather  than  at  the  centroid.  The  ultimate  result  is  that  the 

work  equations  for  model  I  and  model  II  are  identical: 
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3.3  Determination  of  the  Critical  Buckling  Load 


By  setting  the  variation  of  the  work  done  equal  to  the  variation 
of  the  strain  energy  the  value  of  Per  is  once  again  determined.  Since 


the  work  expression  for  the  flat  plate  is  identical  to  that  of  the  I 
beam,  so  too  is  the  equation  for  the  variation  of  the  work: 
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iile  that  for  <5  U  is  slightly  different  owing  to  the  abscence  of  Cw: 
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Equating  the  work  and  strain  energy  and  solving  for  Per  yields: 
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When 


solved  this  quadratic  gives  for  Per  an  expression  which  may  be 


written  in  a  form  similar  to  equation  (2-2)  : 
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(3-1) 


For  typical  stiffeners  the  quantity  ( - — )  will  be  several 

. o  2  6  +  ir“  EI 

orders  of  magnitude  greater  than  ( — — — _.  )  (  —  f—— — Z.  .  Therefore, 

6  +  7T  ^  ^  GJ 

purposes  of  a  design  formula  it  is  sufficiently  accurate  to  ignore 
this  latter  quantity  and  rewrite  equation  (3-1)  as: 
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The  formula  for  Per  for  a  flat  plate  stiffener  subjected  to  a 
point  load  in  the  center  of  the  span  is  given  in  reference  (1)  by 
equation  (6-37),  which  in  this  thesis'  notation  is: 


Per  = 


16.94 
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Now  to  two-decimal  accuracy  equation  (3-2)  may  be  rewritten  as: 

vC 
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The  formulas  are  in  agreement  within  1.3%. 


CHAPTER  4 


MODELS  I  AND  II  -  TORSIONAL  BUCKLIN3  UNDER  COMBINED 
LATERAL  AND  AXIAL  POINT  LOADIN3 

4.1  Description  of  Loading  Condition 

In  chapters  2  and  3  the  failure  by  torsional  buckling  of  a  beam 
subject  to  a  lateral  los  only  has  been  studied;  now  suppose  that  an 
additional  axial  load  is  added  at  the  ends  of  the  beam  as  shown  in 
f igure  (4-1 } . 


Figure  (4-1)  Combined  Lateral  and  Axial  Loading  of  a 


Stiffener 


The  stiffener  is  again  simply  supported  and  laterally  loaded  by  a 
force  P  applied  at  the  upper  surface  in  the  center  of  the  span,  but  an 


additional  force  Q  acts  at  the  centroid  of  each  end  to  compress  the 
stiffener  inward.  The  two  forces  are  acting  at  right  angles  to  each 
other;  P  along  the  Z  axis  and  Q  along  the  X  axis.  This  dual  axis  type 
loading  would  seem  to  be  much  more  representative  of  the  kind  of 
loading  conditions  to  be  encountered  by  the  frames  of  a  ship.  For 
example,  two  commonly  encountered  conditions  by  a  ship  at  sea  in  large 
waves  are  hog  and  sag,  shown  in  figure  (4-2). 
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Here  the  ship  is  supported  by  the  crest  of  a  wave  either  amidships 
(hogging)  or  at  the  ends  (sagging).  In  both  cases  the  vertical 
transverse  frames  of  the  ship,  a  portion  of  which  are  shown  as  dashed 
lines  in  figure  (4-2),  will  be  placed  in  compression,  i.s.  a  Q  force 
applied.  At  the  same  time  then,  the  sides  of  the  ship  will  very 
likely  be  receiving  a  lateral  force  P  from  other  waves  in  a  confused 
sea  or  in  a  worst  case  from  collision  with  another  ship,  floating 
object,  or  possibly  from  within  the  ship  from  unsecured  cargo 
buffeting  about.  How  do  the  two  forces  interact,  and  what  effect  does 
the  addition  of  an  axial  force  have  upon  the  capability  of  the 
stiffener  to  accept  a  lateral  load?  These  are  the  questions  this 
thesis  will  attempt  to  answer  utilizing  the  strain  energy  approach 
which  was  used  in  the  previous  chapters  for  a  lateral  load  alone. 

4.2  Development  of  the  Strain  Energy  Equation 


In  solving  for  the  strain  energy  of  a  stiffener  subjected  to  both 
lateral  and  axial  loads  the  case  of  model  I,  an  I  beam,  will  be 
considered  first;  so  the  energy  equation  will  contain  a  warping 
rigidity  term  Cw: 
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dx 


Again  for  B,  the  angle  of  twist,  the  choice  Bo  cos  I2L  is  made.  The 

d2v  a 

expression  for  - will  however,  no  longer  be  the  same.  Taking  again 

dx2 

an  arbitrary  point  x  as  shown  in  figure  (2-2),  there  will  still  be 
present  to  the  right  of  x  the  force  P/2  creating  a  moment  about  the 
local  £  axis  equal  to  the  expression  in  equation  (2-1): 


There  is  now  though,  another  force  to  the  riqht  of  x  adding  to  this 
moment  about?-  ,  the  axial  force  Q.  Lookinq  down  upon  the  stiffener 
from  above  it  is  seen  that  this  addition  to  M t  is  equal  to  the  force  Q 
times  the  lever  arm  -v  as  shown  in  figure  (4-3): 


Figure (4-3)  Additional  moment  Qv  provided  by 
axial  loading 

The  total  moment  M  t  then  is  now: 

M£  =  |  (|  -  x)  B  +  Qv 


35 


becomes : 


and  the  expression  for 


P  ,  a  v  „  ~x  ,  Qv 
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(4-1) 


We  now  have  a  linear  second  order  differential  equation  which  must  be 
solved  for  v  in  order  that  the  strain  energy  may  be  found.  We  are 
looking  for  a  particular  solution  to  this  differential  equation,  and 
thus  it  will  be  of  the  form  (reference  3): 
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Differentiating  this  expression  for  vn twice  with  respect  to  x  yields: 
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If  v  is  a  solution  to  equation  (4-1)  the  following  is  true: 
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(4-4) 


By  equating  coefficients  of  the  trigonometric  functions  in  equations 
(4-3)  and  (4-4),  four  equations  can  be  found  to  be  solved  for  the  four 


constants  through  .  They  are: 
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Inserting  these  constants  into  equations  (4-2)  and  (4-4)  gives  for  v 

£ 

and  its  second  derivative: 
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For  simplicity  the  coefficients  in  equation  (4-7)  shall  be  labeled 


A,B,  and  C  as  indicated,  and  (4-7)  then  becomes: 
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Now  with  equation  (4-8)  there  is  sufficient  information  to  solve  the 


strain  energy  equation, 
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Carrying  out  the  integration  of  equation  (4-9)  and  reinserting  the 
proper  quantities  for  A,B,  and  C  gives  as  a  final  expression  for 
strain  energy: 
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4.3  Development  of  the  Work  Equation 


The  work  done  will  again  be  found  from  the  principle  that  it  is 
force  times  distance.  The  work  performed  by  the  load  P  is  given  by 
the  same  expression  used  in  chapter  2: 

-  (|  -  x)dx  +  ^  Bo2  (2-3) 

dx^ 

d  2  v 

Now,  however,  there  will  be  a  new  expression  for  — ~  ,  that  of 

dxz 

equation  (4-6) .  Also,  there  must  now  also  be  added  the  work  of  the 


load  Q  at  the  end  x  =  a/2  as  it  moves  through  a  distance  -u  (figure 


Figure  (4-4)  Motion  of  force  Q  through  distance  -u 


The  total  work  done  by  both  loads  will  be: 
a/2 


W  =  P 
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All  the  necessary  quantitites  in  equation 
exception  of  u  _  ,  .  To  find  this  unknown 

strain  (reference  5)  will  be  used: 
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(4-11)  are  known  with  the 
the  definition  of  membrane 


Utilvzing  inextensibility  assumptions  this  membrane  strain  is  set  to 
zero: 
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(4-12) 


As  we  are  dealing  with  failure  of  the  stiffener  by  a  torsional 
(twisting)  mode  augmented  by  an  axial  compressive  force,  it  is 
expected  that  deflection  in  the  Z  direction  will  be  minor  compared 
with  that  in  the  Y  direction.  Therefore,  for  purposes  of  a  design 
relationship  it  will  be  considered  sufficiently  accurate  to  ignore  the 
-4^ —  term  in  equation  (4-12),  and  the  relation  becomes: 
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An  expression  for  u  may  now  be  found  by  solving  equation  (4-13) .  Both 
sides  are  multiplied  by  dx  and  then  integrated: 
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To  find  the  total  work  done  by  Q,  it  is  necessary  to  find  the 
deflection  of  the  extreme  right  hand  side  of  the  beam,  at  x=a/2.  So 
the  limits  of  integration  must  be  over  the  entire  length  of  the  beam, 
x=-a/2  to  x=a/2;  or,  twice  over  half  the  beam,  x=0  to  x=a/2: 
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For  v  equation  (4-5)  is  used,  and  the  first  derivative  of  v  with 
respect  to  x  is: 


41 


4a  Q  +  4EI  tt 
.  z 


(a  Q  EIztt  ) 


2 

a  PBo  *  ttx 

+  (  — 5 - 5-)  cos  — — — 

2a  Q  +  2  El  ii 
v  z 


,  aTTPBo  x  _  ttx 

(  — 5 - 5-)  x  sin  — — 

2a^Q  +  2  El  tt  a 

2 


Again  for  simplicity  the  coefficients  in  equation  (4-15)  are  labeled 
D,E,F,  and  G  as  shown;  and  equation  (4-14)  becomes: 
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Carrying  out  the  integration  of  equation  (4-16)  and  replacing 


coefficients  gives  for  the  displacement  U  x=a/2: 
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All  quantities  are  now  known  to  determine  the  work  done  by  the 


combined  lateral  and  axial  loads.  Utilyzing  the  proper  expressions 
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for  B,  (eq.  4-6),  and  u  (eq.  4-17),  equation  (4-11)  becomes: 
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Integrating  equation  (4-18)  yields  for  work: 


W  = 


2  2  3 

P  Bo  a 


El 


2  2 

,  6  +  TT  ,  ,  6  +TT  \  , 

(  - — )  -  (  —  T  '  ( 


a20 


96  it 


24tt 


2  2 
4a  Q  +  4EIztt 


+  ( 


2  2  5 

P  ho  a  Q 


(4a2Q  +  4EIztt2)2 


2  , 

TT  -  6 


12 


)  -  ( 


24EI  it 
z 


4a2Q  +  4EI  tt2 
z 


8EI  V 
z 


4a2Q  +  4EI  tt  2 


Ph  n.2 
-: — 30 


+  ( 


(4-19) 


By  conservation  of  energy  6W  =  <5U;  thus  equations  (4-20)  and  (4-21) 
may  be  equated,  and  the  ultimate  relation  between  P  and  Q  appears  as 
the  following  quadratic  with  respect  to  P: 
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Equation  (4-22)  presents  the  relationship  between  P  and  Q  when  they 
are  applied  to  an  I  beam.  If  both  fprces  are  applied  in  like  manner 
to  model  II,  the  flat  bar  stiffener,  the  only  difference  in  equation 

IT4 

(4-22)  is  that  the  final  term,  ECw ( - ,  will  be  absent  since  Cw 

4aJ 


for  the  flat  bar  is  zero. 


CHAPTER  5 


EVALUATION  OP  RESULTS  OF  THE 
BUCKLINS  LOAD  EQUATION 


5.1  Nondimensionalization  of  the  Buckling  Load  Equation 

It  is  now  our  desire  to  use  equation  (4-22)  to  study  the 

i 

interaction  between  the  lateral  and  axial  loads  P  and  Q  and  determine 
how  changes  in  the  physical  dimensions  of  a  stiffener  effect  its 
ability  to  withstand  against  failure.  The  most  interesting  comparison 
to  be  made  from  equation  (4-22)  would  be  a  plot  of  lateral  load  P 
versus  axial  load  Q,  i.e.  a  means  to  predict  for  a  given  amount  of 
lateral  load  how  much  simultaneous  axial  load  a  stiffener  will  take 
before  buckling  and  vice  versa. 

In  order  to  construct  such  a  plot  it  will  be  convenient  to  first 
nond imensionalize  equation  (4-22)  .  An  ideal  means  to 
nondimensionalize  P  and  Q  is  to  express  them  as  a  fraction  of  the 
Euler  buckling  load  (equation  2-5  of  reference  1) : 


a 
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The  Euler  critical  buckling  load  Qe  is  that  load  applied  as  in  the 
manner  shown  in  figure  5-1  that  will  just  cause  in-plane  deflection  of 
the  beam  without  the  twisting  effects  of  torsion. 


Figure  (5-1)  Euler  Buckling  of  a  Beam 


Use  of  the  Euler  buckling  load  as  a  nondimensionaliz ing  parameter  i 
particularly  well  suited  to  equation  (4-22)  noticinq  the  frequent 
occurrence  in  it  of  the  expression 

4a2Q  +  4EI  y2 


which  may  be  written,  if  0  is  taken  as  positive  in  tension  as: 

-4a2Q  +  4a2Qe  =  4a2Qe  (  1  -  ) 

Qe 


The  final  dimensionless  form  of  (4-22)  will  be: 
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where 


L  =  1  - 


To  three  decimal  accuracy  equation  (5-1)  becomes: 


'of)2  f  -816  -■§?  [  i-*33  <E 
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)  2  +  3.701  (-r)  3  -  2.467  (})  4 

Li  L 


949  ^  “  24.352  ( — ^2 — )  =  0 

z  a  l 
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+  2-467(^1  | 


(5-2) 


Equation  (5-2)  will  be  applicable  to  an  I  beam,  and  to  make  it  useful 
for  a  flat  bar  stiffener  as  usual  set  Cw  to  zero  and  delete  the  final 
term. 


5.2  Model  I  -  Effects  of  Varying  Stiffener  Dimensions  on 
the  Value  of  the  Critical  Buckling  Load 


The  results  of  equation  (5-2)  for  I  beams  of  various  sizes  are 
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shown  graphically  in  figures  (5-2)  to  (5-5).  The  points  on  the  graphs 
were  derived  from  computer  analysis  shown  in  Appendix  A.  Along  the 
abscissa  of  each  of  the  graphs  is  shown  the  dimensionless  ratio  of 
axial  load  Q  to  the  Euler  buckling  load  Qe,  while  the  ordinate  gives 
the  ratio  of  the  critical  torsional  buckling  load  Per  to  the  Euler 
load.  Thus  by  calculating  the  Euler  buckling  load  of  a  stiffener  and 
being  given  a  value  of  axial  load  Q,  we  can  enter  the  graphs  with  a 
value  of  Q/Qe  and  obtain  a  prediction  of  how  much  lateral  force  Per, 
applied  at  the  center  of  the  span,  that  the  stiffener  will  withstand 
before  torsional  buckling  occurs.  For  example,  in  figure  (5-2)  if  an 
axial  force  equal  to  a  value  of  .4Qe  is  applied  to  stiffener  3  with 
h/a  ratio  equal  to  .14,  then  a  lateral  force  of  just  .16  Qe  should 
produce  failure  of  the  stiffener  by  torsional  buckling. 

Though  no  previous  design  formulas  relating  simultaneous  axial 
and  lateral  loading  of  a  stiffener  to  cause  torsional  buckling  were 
found,  the  results  of  equation  (5-2)  do  appear  reasonable.  It  is 
possible  to  verify  the  predicted  results  of  Per /Qe  at  the  intersection 
with  the  ordinate  (Q/Qe=0).  For  example,  again  for  stiffener  3  of 
figure  (5-2)  with  h/a  equal  to  .14,  the  complete  set  of  scantlings  for 
this  I  beam  are: 

a*3 .07m 
h* . 43m 
b= .09m 
t  f = . 02m 
tw= . 007m 
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Numbers  adjacent  to  each 
curve  correspond  to  stiffener 
numbers  in  Appendix  A 


h/a  increasing 


Figure  (5-2)  Pcr/Qe  vs  O/Qe  for  I  beam  varying  h/a 


fi! 


Utilyzing  these  dimensions  the  warping  constant,  St.  Venant's  torsion 

constant,  and  moment  of  inertia  of  the  stiffener  about  the  web  are: 

c„  .  (.02m)  ( .  4  3m)  ^  ( .  09m)  ^  .  la2J  x  10*7  m6 

24 
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u  —  * 5  *  2  9  2,  x  1 0  rn 


I  =  ■ Sisl  Li02] m]  _  2.430  X  10_6m4 

z  12 


Now  the  Euler  buckling  load  Qe  for  this  beam  calculated  from  equation 

(2-5)  of  reference  (1)  is: 
it2EI 

Qe  =  - 

a 

tt2  (2.068  x  108  KN/m2)  (2.430  x  10~6m4) 

(3.07m) 2 


=  526.2  KN 


To  find  the  critical  load  for  torsional  buckling  from  a  lateral  load 
alone  we  use  equation  (6-18)  of  reference  (1) : 


Per  = 


If  is  obtained  by  entering  figure  (2-4)  of  this  thesis  with 

(3.07m)2 (7.954  x  IQ7  KN/m2) (5.292  x  10~7m4) 
(2.068  x  108  KN/m2) (1.123  x  10_7m6) 


2  GJ 
a  ECw 


and  reading  off  of  15.3.  Now  the  critical  torsional  buckling  load 


=  236.1  KN 

A  calculated  ratio  of  Per /Qe  then  is: 

Per  _  236.1  KN 
Qe  526.2  KN 

=  .449 

The  value  of  Pcr/Qe  obtained  from  equation  (5-2)  for  this  same 

stiffener  and  shown  on  figure  (5-2)  for  Q/Qe=0  is  .455;  an  agreement 

of  1.3%.  Table  5-1  lists  the  values  for  Pcr/Qe  at  Q/Qe=0  for  all  the 

curves  of  figures  (5-2)  through  (5-5)  and  shows  the  good  correlation 

between  the  results  of  equation  (5-2)  and  the  values  calculated  using 

equations  (2-5)  and  (6-18)  of  reference  (1).  All  cases  are  within  2% 

agreement  except  stiffener  1,  figure  (5-2),  and  stiffener  9,  figure 

(5-4).  In  these  two  cases  a  problem  arises  in  trying  to  pick  off  an 

accurate  value  of  ^  from  figure  (2-4)  for  use  in  equation  (6-18) .  We 

are  in  the  vicinity  of  a  GJ  equal  to  one  on  figure  (2-4)  .  Here  the 

ECw 

values  of  If  are  changing  quite  rapidly,  and  a  small  error  in  It  can 
lead  to  an  error  of  several  percentage  points  in  Pcr/Qe. 


Curve 


Pcr/Qe  at  Q/Qe=0 
equation  (5-2) 


Pcr/Qe  AT  Q/Qe=0 
reference  (1) 


Percent 


Agreement 


Fig 

5-2  ,Beam 

1 

.986 

1.056 

6.6% 

Beam 

2 

.586 

.583 

.5% 

Beam 

3 

.455 

.449 

1.3% 

Fig 

5-3 , Beam 

4 

.596 

.598 

.3% 

Beam 

5 

.586 

.583 

.5% 

Beam 

6 

.573 

.575 

.5% 

Fig 

5-4  ,Beam 

7 

.945 

.938 

.7% 

Beam 

8 

.586 

.583 

.5% 

Beam 

9 

.486 

.515 

5.6% 

Fig 

5-5 ,Beam 

10 

.748 

.734 

1.9% 

Beam 

11 

.586 

.583 

.5% 

Beam 

12 

.575 

.587 

2.0% 

Table  (5 

-1) 

Comparison  of 

values  for  Pcr/Qe 

at  Q/Qe=0  for  an  I 

beam 

i  as  calculated 

by  equation 

(5-2)  versus  values 

calculated  by 

equations  (2-5)  and  (6-18)  of  reference  (1). 


At  the  right  hand  side  of  the  curves  of  figures  (5-2)  through 
(5-5)  the  graphs  asymptotically  approach  the  x  axis,  but  have  no  value 
at  Q/Qe  =  1.  This  is  because  at  Q/Qe=l  the  expression  1/L  in  equation 
(5-2)  will  have  the  indeterminate  form  1/0.  This  is  to  be  expected 


though,  because  at  Q/Qe=l  the  stiffener  has  failed  due  to  Euler 
buckling,  and  equation  (5-2)  can  no  longer  be  used  to  describe  the 
relation  between  Per  and  Q. 

The  general  trend  of  all  the  curves  is  similar,  a  rapidly 
decreasing  initial  slope  with  the  values  for  Pcr/Qe  falling  off 
approximately  50%  by  Q/Qe  =.37,  then  tapering  off  of  the  slope  as  the 
curves  approach  Q/Qe=l.  The  dimensionless  parameters  that  were  varied 
in  constructing  the  graphs  were  h/a,  h/tw,  b/tf,  and  tf/tw.  In  each 
figure  the  varied  parameter  has  been  doubled  twice  to  construct  the 
three  curves.  It  is  seen  that  the  ratios  which-  exhibit  the  greatest 
influence  on  Pcr/Qe  are  h/a  and  b/tf.  Varying  h/tw  has  practically  no 
effect  upon  Pcr/Qe,  while  changing  tf/tw  at  first  appears  to  have  a 
fairly  significant  effect  but  seems  to  steady  out  after  tf/tw  gets 
above  about  2.8. 

5.3  Model  II  -  Effects  of  Varying  Stiffener  Dimensions  on  the 

Value  of  the  Critical  Buckling  Load 

Figure  (5-6)  shows  the  results  of  equation  (5-2)  for  model  II, 
the  flat  bar  stiffener.  Once  again  the  values  of  Pcr/Qe  for  Q/Qe=0 
may  be  checked  against  those  obtained  using  the  formulas  of 
Timoshenko.  Consider  for  example  stiffener  13  with  h/a=.0523.  The 
stiffener  will  have  dimensions  and  factors  as  follows: 
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a  =  3.059m 


h  =  .16m 


tw  =  .0064m 


J 


I 

z 


( . 16m)  ( . 0064m) 9 
3 

(.16m) (,0064m)3 
12 


Cw  =  0 

=  1.398  x  10-8m4 

3.495  x  10~9  m4 


The  Euler  buckling  load  from  equation  (2-5)  of  reference  (1)  is: 


Qe  = 


tt2  (2.068  x  108KN/m2)  (3.495  x  10~9m4) 
(3.059m) 2 


=  . 762KN 


Per  for  the  flat  bar  from  equation  (6-37)  of  reference  (1)  is: 
16.94  ^ El  GJ 


Per  = 


(1 


-  *87  5  / 


El 


GJ 


) 


=1.564  KN 


Per  _  1.564  KN 

Qe  .762  KN 


2.052 


The  value  of  Pcr/Qe  obtained  from  equation  (5-2)  is  2.079,  an 
agreement  of  1.3%.  The  agreement  of  the  values  for  the  other  two 
curves  are  shown  in  table  (5-2)  .  It  is  seen  that  there  is  good 


agreement  in  all  three  cases. 


Curve  Pcr/Qe  AT  Q/Qe=0  Per /Qe  AT  Q/Qe=0  Percent 

1 

equation  (5-2)  reference  (1)  Agreement 

FIG  5-6, Beam  13  2.079  2.051  1.3% 

FIG  5-6, Beam  14  2.005  1.978  1.4% 

FIG  5-6, Beam  15  1.864  1.817  2.6% 

TABLE  (5-2)  Comparison  of  values  of  Pcr/Qe  at  Q/Qe=0  for  a  flat  bar 
as  calculated  by  equation  (5-2)  versus  values  calculated  by  equation 
(6-37)  of  reference  (1). 


As  before  with  the  I  beam,  equation  (5-2)  is  undefined  for  the 
flat  bar  at  Q/Qe=l  owing  to  the  quantity  1/L  being  equal  to  1/0;  but 
this  is  acceptable  since  the  bar  has  undergone  Euler  buckling  at 
Q/Qe=l.  The  shape  and  trends  of  the  curves  are  identical  to  the  I 
beam,  and  again  the  curve  falls  50%  by  about  Q/Qe=.36.  It  is 
interesting  to  note  that  whereas  for  the  I  beam  Pcr/Qe  is  always  less 
than  one,  Pcr/Qe  for  the  flat  plate  starts  out  greater  than  one  by  a 
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L6:< 


factor  of  approximately  two  and  remains  greater  than  one  until 
Q/Qe=.29.  This  would  indicate  the  much  lower  relative  Qe  for  the  flat 
bar  and  its  much  greater  susceptibility  to  failure  by  Euler  buckling 
than  the  I  beam. 

Owing  to  the  geometric  simplicity  of  the  flat  bar,  the  only 
dimensionless  parameter  that  may  be  varied  independently  is  h/a.  If 
h/tw  is  varied  while  keeping  h/a  constant,  there  will  be  no  change  in 
equation  (5-2) .  As  Cw  is  zero  and  the  ratio  J/Iz  for  a  flat  bar  is 
always  4,  the  only  term  that  will  alter  equation  (5-2)  is  h/a. 
Moreover,  as  can  be  seen  from  figure  (5-6),  even  the  changing  of  h/a 
has  very  little  effect  upon  Pcr/Qe. 
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CHAPTER  6 


CONCLUSIONS 

In  this  thesis  has  been  studied  the  problem  of  torsional  buckling 
of  I  beam  and  flat  bar  type  stiffeners  under  combined  lateral  and 
axial  loads.  First  examined  was  the  case  of  lateral  loading  alone. 

An  energy  method  was  employed  to  determine  critical  buckling  stress 
Per  of  a  stiffener,  and  equations  (2-4)  and  (3-1)  were  developed  to 
cover  lateral  loading  of  I  beams  and  flat  bar  stiffeners  respectively. 
Good  agreement  was  shown  between  equations  (2-4)  and  (3-1)  and  the 
published  formulas  of  Timoshenko.  One  distinct  advantage  of  equation 
(2-4)  for  an  I  beam  is  that  it  gives  an  exact  formula  for  the 
coefficient  Y  in  reference  (1)  equation  (6-18)  rather  than  a  tabular 
listing  of  Y  for  different  sizes  of  stiffeners.  As  was  noted  in 
chapter  5,  when  trying  to  verify  some  of  the  values  of  Pcr/Qe  it 
becomes  qux^e  difficult  to  pick  off  an  accurate  value  of  Y  when  near 
the  bottom  of  the  curve  in  figure  (2-4).  In  particular,  for  the  two 
cases  in  Table  (5-1)  where  the  percent  agreement  between  Pcr/Qe  as 
calculated  from  this  thesis'  equations  and  those  of  reference  (1)  were 
6.6%  and  5.6%,  the  differences  fall  to  just  .5%  and  .2%  when  Y  is 
calculated  from  equation  (2-4)  rather  than  read  from  the  graph. 


62 


Next  was  approached  the  problem  of  combined  lateral  and  axial 
loading  of  a  stiffener.  Again  an  energy  approach  to  the  problem  was 
used,  and  the  ultimate  result  was  equation  (4-22)  which  when 
nondimensionalized  using  Euler  buckling  load  took  the  form  of  equation 
(5-1).  Though  no  other  previous  formula  relating  torsional  buckling 
under  combined  loads  was  found  for  comparison,  the  end  points  of  the 
curves  of  figures  (5-2)  through  (5-6)  could  be  checked  and  were  found 
to  be  in  good  agreement  with  results  obtained  using  other  published 
formulas.  The  trends  of  all  the  curves  appear  reasonable:  a  constant 
decrease,  more  rapid  at  first,  of  the  amount  of  lateral  load  a 
stiffener  can  sustain  as  an  increasing  axial  load  is  applied.  As 
another  vote  of  confidence,  if  0  is  set  to  zero  in  equation  (4-22)  the 
expression  reverts  to  that  of  equation  (2-4)  which  was  verified 
against  equation  (6-18)  of  reference  (1) .  An  area  for  research  would 
be  to  experimentally  verify  the  results  of  equation  (4-22),  though  it 
would  certainly  be  expensive  due  to  the  complexity  of  a  test  apparatus 
capable  of  imposing  and  measuring  simultaneous  lateral  and  axial 
loads. 

Equation  (4-22)  is  somewhat  lengthy  for  hand  calculations, 
particularly  when  the  constants  Cw,  H,  and  Iz  must  also  be  found 
before  using  it.  It  is  not  too  cumbersome  if  only  one  stiffener  is 
involved,  but  if  a  number  of  stiffeners  are  being  investigated  it  does 
lend  itself  well  to  a  computerized  study,  especially  if  an 
investigation  is  to  be  done  varying  stiffener  dimensions.  The  simple 
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BASIC  program  with  accompanying  results  shown  in  Appendix  A  was 
written  and  utilyzed  to  generate  the  values  for  constructing  figure 
5-2  through  5-6. 
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APPENDIX  A 


Computer  Adaptation  of  Equation  (4-22) 
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